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ABSTRACT 

Wc present detailed evolutionary calculations for carbon - oxygen - and helium - core 
white dwarf models with masses ranging from M = 0.1 to M = 1.2 M® and for 
metallicities Z — 0.001 and Z = 0. The sequences cover a wide range of hydrogen en- 
velopes as well. We have taken finite - temperature effects fully into account by means 
of a detailed white dwarf evolutionary code, in which updated radiative opacities and 
equations of state for hydrogen and helium plasmas are considered. The energy trans- 
port by convection is treated within the formalism of the full - spectrum turbulence 
theory, as given by the self - consistent model of Canuto, Goldman and Mazzitclli. 
Convective mixing, crystallization, hydrogen burning and neutrino energy losses are 
taken into account as well. 

The set of models presented here is very detailed and should be valuable, par- 
ticularly for the interpretation of the observational data on low - mass white dwarfs 
recently discovered in numerous binary configurations and also for the general problem 
of determining the theoretical luminosity function for white dwarfs. In this context, 
we compare our cooling sequences with the observed white dwarf luminosity function 
recently improved by Leggett, Ruiz & Bergeron and we obtain an age for the Galactic 
disc of ss 8 Gyr. Finally, we applied the results of this paper to derive stellar masses 
of a sample of low - mass white dwarfs. 

Key words: evolution - stars: stars: interiors - luminosity function, mass function - 
pulsars: general - white dwarfs 



1 INTRODUCTION 

Numerous observations carried out over recent years have 
presented strong evidence that low - mass, helium white 
dwarf stars are the product of the evolution of certain close 
binary systems. Indeed, low - mass white dwarfs have been 
detected in binary systems containing, for instance, another 
white dwarf (Marsh 1995; Marsh, Dhillon & Duck 1995; 
Marsh & Duck 1996; Moran, Marsh & Bragaglia 1997), 
millisecond pulsar (Lundgren et al. 1996; see also Backer 
1998) or a yellow giant (Landsman et al. 1997). In particu- 
lar, Moran et al. (1997) found the binary system WD 0957 
- 666 (consisting of two low - mass white dwarfs) to have 
an orbital period of only 1.46 h, which is short enough for 
the binary to merge within only 2.0 x 10 8 yr. Very recently, 
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Edmonds et al. (1998) have reported the presence of a candi- 
date helium white dwarf in the globular cluster NGC 6397. 
On theoretical grounds, recent population models of close 
binaries (Iben, Tutukov & Yungelson 1997) suggest a high 
probability of discovering helium white dwarfs in close bina- 



Detailed evolutionary models of low - mass white dwarfs 
may provide valuable information not only on the white 
dwarf itself but also on the companion object and even on 
the past evolution of the system (see, for instance, Burderi, 
King & Wynn 1996 and Hansen & Phinney 1998a). In this 
regard, the analysis carried out, notably by van Kerkwijk, 
Bergeron & Kulkarni (1996), is worth mentioning. Indeed, 
from spectroscopic data inferred from its low - mass white 
dwarf companion combined with a theoretical mass - ra- 
dius relation for the white dwarf, these authors found the 
mass of the pulsar PSR J1012 + 5307 to be between 1.5 and 
3.2 Mq . Needless to say, detailed models of helium white 
dwarfs, together with further observations, are needed in 
order to achieve a more precise determination of the pul- 
sar mass and hence to constrain the equation of state at 
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the high densities appropriate for neutron stars. In addition, 
an independent determination of the age of many millisec- 
ond pulsars can be inferred from the study of the cooling 
of their helium white dwarf companions, which is valuable 
for understanding the nature and origen of such systems. 
Another strong motivation for constructing improved white 
dwarf evolutionary sequences is the fact that, thanks to the 
Hubble Space Telescope, it has been possible to detect the 
low - luminosity tail of the white dwarf population in glob- 
ular clusters. Accordingly, white dwarf evolutionary tracks 
would provide an independent way of measuring the age and 
distance of such clusters (see, e.g., Richer et al. 1995; Von 
Hippel, Gilmore & Jones 1995 and Renzini et al. 1996). 

In view of these considerations, we present in this pa- 
per new grids of white dwarf evolutionary models for dif- 
ferent hydrogen envelopes and stellar masses. The emphasis 
is placed mainly on low - mass, helium white dwarfs, the 
detailed study of which has recently began to be under- 
taken. As a matter of fact, Althaus & Benvenuto (1997a) 
and Benvenuto & Althaus (1998) carried out an analysis 
of the structure and evolution of low - mass white dwarfs 
based on a updated physical description, such as new opac- 
ities and equations of state, and the employment of a new 
convection model more physically sound than the mixing 
length theory. In a still more recent study, Hansen & Phin- 
ney (1998a) presented evolutionary calculations for these ob- 
jects as well. However, the evolutionary sequences for their 
more massive models do not converge to the Hamada & 
Salpeter (1961) predictions for zero - temperature, pure - 
helium configurations, thus resulting in models with under- 
estimated surface gravities. This can be seen from figure 16 
of Hansen & Phinney (1998a). Note that the surface gravity 
for their more massive models with a hydrogen envelope of 
Me/Mq = 1CF 6 is substantially lower than the Hamada - 
Salpeter values. Such a discrepancy cannot be attributed to 
the hydrogen layer, since a thin hydrogen envelope intro- 
duces a very small correction to the stellar radius of zero - 
temperature, pure helium models (see Benvenuto & Althaus 
1998). 

By contrast, the study of evolution of carbon - oxy- 
gen white dwarfs have captured the interest of numerous 
investigators such as Lamb & Van Horn (1975), Iben & Tu- 
tukov (1984), Koester & Schonberner (1986), D'Antona & 
Mazzitelli (1989), Tassoul, Fontaine & Winget (1990), Wood 
(1992), Benvenuto & Althaus (1997) and Althaus & Ben- 
venuto (1998). In particular, Iben & Tutukov (1984) were 
the first in computing evolutionary models of white dwarfs 
with hydrogen burning, showing that hydrogen burning in 
cooling white dwarfs could be an important energy source. 

With the calculations that we present here we amply ex- 
tend those presented in Althaus & Benvenuto (1997a) and 
Benvenuto & Althaus (1998), in which the effects of convec- 
tion, neutrino losses and different hydrogen envelopes on the 
structure and evolution of helium white dwarfs were care- 
fully analysed (we should mention that the results shown 
in Benvenuto & Althaus 1998 correspond to a metallicity 
of Z « and not to Z — 0.001, as stated in that work). 
Furthermore, we extend our calculations to the case of more 
massive carbon - oxygen white dwarfs. Our grid for carbon 
- oxygen models is likewise very detailed, which may be of 
relevance in the study of, for instance, the general prob- 
lem of determining the theoretical white dwarf luminosity 



function and the assessment of the age of the Galactic disc. 
This subject has been recently addressed by Leggett, Ruiz 
& Bergeron (1998), who have greatly improved the deter- 
mination of the observed luminosity function for cool white 
dwarfs. In this regard, we shall derive theoretical luminosity 
functions from our cooling sequences in order to compare 
with the Leggett et al. observational data. 

The results presented here constitute a very detailed 
and updated set which will be suitable, for instance, for the 
interpretation of recent and forthcoming observational data 
about low - mass white dwarfs in close binary systems. Fi- 
nally, we applied our evolutionary models with helium cores 
to derive stellar masses of a sample of low - mass white 
dwarfs. 



2 COMPUTATIONAL DETAILS 

The evolutionary sequences have been obtained with the 
same evolutionary code and input physics we employed in 
our previous works on white dwarf evolution, and we refer 
the reader to Althaus & Benvenuto (1997a, 1998) and Ben- 
venuto & Althaus (1998) as well as to the references cited 
therein for details. In what follows we restrict ourselves to a 
few brief comments. 

The code has been written following the method pre- 
sented by Kippenhahn, Weigert & Hofmeister (1967) for 
calculating stellar evolution. In particular, to specify the 
surface boundary conditions we perform three envelope in- 
tegrations (at constant luminosity) from photospheric start- 
ing values inward to a fitting mass fraction M\/M ~ 10 -16 , 
where Mi corresponds to the first Henyey mass shell and M 
is the total mass of the white dwarf model. In our code the 
value of Mi is automatically changed over the evolution so 
as to keep the thickness of the envelope as small as possible. 
This provides an accurate description of the outer layers of 
our white dwarf models. The interior integration is treated 
according to the standard Henyey technique as described by 
Kippenhahn et al. (1967). 

The constitutive physics of our code is as detailed and 
updated as possible. Briefly, for the low - density regime, we 
consider the equation of state of Saumon, Chabrier & Van 
Horn (1995) for hydrogen and helium plasmas. The treat- 
ment for the high - density, completely ionized regime ap- 
propriate for the white dwarf interior is based on our own 
equation of state. This includes ionic and photon contribu- 
tions, coulomb interactions, partially degenerate electrons, 
quantum corrections for the ions and electron exchange and 
Thomas - Fermi contributions at finite temperature (see Al- 
thaus & Benvenuto 1997a for details). Radiative opacitites 
for the high - temperature regime (T > 6000 K) are those 
of OPAL (Iglesias & Rogers 1993), whilst for lower temper- 
atures we use the Alexander & Ferguson (1994) molecular 
opacities (or the Cox & Stewart 1970 tabulation for pure he- 
lium composition). Two extreme values for metallicity have 
been considered in the envelope: Z= and Z= 0.001. We 
should mention that, owing to the lack of reliable low - tem- 
perature (T < 6000-ftT) opacities for helium composition, our 
low - luminosity models with helium atmospheres should be 
regarded with caution, particularly their ages. Conductive 
opacities for the liquid and crystalline phases and the vari- 
ous mechanisms of neutrino emission relevant to white dwarf 
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Table 1. Available white dwarf evolutionary sequences with hydrogen envelopes 
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This table shows available evolutionary sequences for white dwarf models with a hydrogen envelope. We 
use an abbreviated notation to indicate the core composition, the stellar mass in tenths of solar mass units, 
the mass fraction of the hydrogen envelope and the envelope metallicity. For instance, He252e3 .z0 stands 
for an evolutionary sequence of 0.25Mq models with a helium core composition, a hydrogen envelope of 
Mjj/M = 2 X 10~ and an envelope metallicity of Z = 0. OC means a oxygen - carbon core composition 
(Note that OC sequences for model more massive than IMq are indicated with the same notation). We 
also provide a column (T* ff ) with the effective temperature (in thousand K degrees) at which each sequence 
starts. 



interiors are taken from the works of Itoh and collaborators 
(see Althaus & Benvenuto 1997a for details). We also include 
in our code the complete network of thermonuclear reaction 
rates for hydrogen burning corresponding to the proton - 
proton chain and the CNO bi - cycle. Nuclear reaction rates 
are taken from Caughlan & Fowler (1988) and j3 - decay 
rates from Wagoner (1969). Electron screening is from Wal- 
lace, Woosley & Weaver (1982). We use an implicit method 
of integration to compute the change of the following chem- 
ical species: 1 H, 2 H, 3 He, 4 He, 7 Li, 7 Be, 8 B, 12 C, 13 C, 13 N, 
14 N, 15 N, 15 0, 16 0, 17 and 17 F. 

Another important feature of our evolutionary se- 
quences is that the energy transport by convection is de- 
scribed by the full - spectrum turbulence theory (see Canuto 
& Mazzitelli 1991, 1992 and references cited therein for de- 
tails) , which represents a great improvement compared with 
the mixing length theory of convection used thus far in most 
of white dwarf studies. As a matter of fact, the Canuto & 



Mazzitelli theory, which has successfully passed a wide va- 
riety of laboratory and astrophysical tests (Canuto 1996), 
takes into account the whole spectrum of turbulent eddies 
necessary to compute the convective flux accurately in the 
almost inviscid stellar interiors. For the set of sequences pre- 
sented in this paper we have considered the recent improve- 
ment to this convection theory introduced by Canuto, Gold- 
man & Mazzitelli (1996), which has been shown to provide a 
good agreement with recent observational data on pulsating 
white dwarfs (Althaus & Benvenuto 1997b) . The model pre- 
sented by Canuto et al. improves over its predecesor (where 
the rate of input energy is given by the linear growth rate) 
in the fact that the growth rate is computed as a function of 
the turbulence itself, thus ensuring a self - consistent treat- 
ment. At intermediate and low convective efficiency , this 
feature leads to larger convective fluxes as compared with 
the Canuto & Mazzitelli (1992) model. It is worthwhile to 
mention that the mass - radius relation and ages correspond- 
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Figure 1. The location of the top and the base of the convec- 
tion zone in terms of the outer mass fraction q (q = 1 — M r /M) 
versus T e ff for a 0.3 Mq white dwarf model with a hydrogen 
envelope of Mr/M = 10 — 3 according to different theories of con- 
vection and metallicities. It is clear that the mixing temperature 
for models with thin hydrogen envelopes will be dependent upon 
the assumed theory of convection. At low T e f t values, the depth 
reached by the base of the convection zone is independent of the 
treatment of convection. Note also the deeper final extent reached 
by convection in the case of metallicity Z = 0. 

ing to our white dwarf models are practically insensitive to 
the convection theory employed. In contrast, the size of the 
outer convection zone in an intermediate effective tempera- 
ture (T e ff ), evolving white dwarf is strongly dependent upon 
the assumed treatment of convection. Hence, a trustworthy 
model of stellar convection is to be employed to get reliable 
T c ff values at which thin hydrogen envelopes mix with the 
underlying helium (see Benvenuto & Althaus 1998). In this 
context, our calculations represent an improvement over pre- 
vious white dwarf studies based on the mixing length theory 
of convection. To clarify this point better, we show in Fig- 
ure ^ the behaviour of the evolving outer convection zone in 
terms of T c g for our 0.3 Mq model with a thick hydrogen 
envelope. In addition to Canuto et al.'s results, we include in 
the figure the predictions given by the ML1, ML2 and ML3 
versions of the mixing length theory amply used in white 
dwarf studies (see, e.g., Tassoul et al. 1990). It is clear that 
the mixing T e « for models with thin hydrogen envelopes 
depends upon the convection theory. Another observation 
that we can make from this figure is that the base of the 
convection zone (for the case Z — 0) ultimately reaches a 
final extent at an outer mass fraction of 2 x 10 -4 (6 x 10 -5 
Mq ), irrespective of the treatment of convection. This re- 
sult is in good agreement with the predictions of Hansen & 
Phinney (1998a) for the same model. Note also the larger 
final extent of convection zone for lower metallicity, which 
as we shall see, gives rise to considerable differences in the 
evolutionary times at low luminosities. 
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Figure 2. Chemical profile for our carbon - oxygen core models 
versus the fractional mass. Solid lines correspond to carbon and 
medium dashed lines to oxygen. 



Our white dwarf initial models of different masses and 
hydrogen envelopes have been obtained following the artifi- 
cial evolutionary procedure described by Benvenuto & Al- 
thaus (1998). The carbon - oxygen core models all have the 
same core chemical composition profile shown in Figure ^| 
This chemical profile was calculated by D'Antona & Mazz- 
itelli (1989) for the progenitor evolution of a 0.55 Mq white 
dwarf. We adopt this profile for all of our models, in spite 
of the changes that are expected to occur for more massive 
models as a result of differences in the evolution in progeni- 
tor objects. We would need, in order to improve this assump- 
tion, detailed calculations of the pre - white dwarf evoltution 
of these objects. To our knowledge, such calculations are not 
available. Because of the fact that the mass of the hydrogen 
envelope in white dwarfs is poorly constrained by theoretical 
calculations of the pre - evolution of these objects, particu- 
larly in the case of helium white dwarfs where the uncertain- 
ties regarding the mass exchange episodes are more severe, 
we decide to treat the mass of the hydrogen envelope as es- 
sentially a free parameter. It is worthwhile to mention that 
our evolving low - mass white dwarf models should be con- 
sidered as evolutionary stages that can be asymptotically 
reached by helium white dwarfs resulting from close binary 
evolution. In this study we have not computed such binary 
evolution, and we refer the reader to Iben & Livio (1993) 
for a review. Needless to say, the starter model choice af- 
fects the initial evolution of all of our models, particularly 
the age (see Althaus & Benvenuto 1997a for details). 

In closing, we have included in our calculations the re- 
lease of latent heat during crystallization (see Benvenuto & 
Althaus 1997 for details) and convective mixing. 
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Table 2. Available white dwarf evolutionary sequences without a hydrogen envelope 



Sequence 


1 B B 


Sequence 


n i i 
J cff 


Sequence 


J eff 


Sequence 


J eff 


neuyyu.zo 


P. fi9 


neiuoy .zo 


/ .uu 


nGlUy 1 .Zo 


7 AX. 


u i i Aa „q 

xiei i^o.zo 


7 Q 1 ^ 


XlClZUO.ZO 


£ A*\ 


ITr>1 9f^ -7^ 

xieizoo.zo 


q no 
y.uz 


xieiozo.zo 


Q £9 


neioyo.zo 


1 n 9R 

1U. ZD 


XlCl^iDZ.ZO 


i n qq 
ru.yy 


xieiouu.zo 


1 1 


U i c r 7t; „q 


19 17 


XieiOO'ri-ZO 


1 9 on 
iz.yu 


HclYob.zo 


13.51 


Helozo.zo 


14.08 


TT — 1 m ,1 1 

Hel914.zo 


14.82 


riczUlU.zo 


15.85 


He2110.z3 


17.26 


He2216.z3 


18.27 


He2327.z3 


20.01 


He2443.z3 


21.42 


He2565.z3 


22.83 


He2693.z3 


24.61 


He2829.z3 


25.29 


He2970.z3 


26.70 


He3118.z3 


26.44 


He3274.z3 


26.77 


He3438.z3 


27.90 


He3610.z3 


28.07 


He3790.z3 


30.09 


He3979.z3 


29.63 


He4179.z3 


30.60 


He4388.z3 


32.49 


He4607.z3 


34.63 


He4837.z3 


37.30 


He5079.z3 


37.69 






OC5000.z3 


57.61 


OC5200.z3 


47.60 


OC5400.z3 


51.96 


OC5600.z3 


54.21 


OC5800.z3 


59.70 


OC6000.z3 


78.72 


OC7000.z3 


78.59 


OC8000.z3 


93.85 


OC9000.z3 


94.82 


OC1000.z3 


101.57 


OC1100.z3 


107.97 


OC1200.z3 


100.01 



This table shows available evolutionary sequences for white dwarf models without a hydrogen envelope. 
We use an abbreviated notation to indicate the core composition, the stellar mass in tenths of solar mass 
units and the envelope metallicity. For instance, He3118.z3 stands for an evolutionary sequence of 0.3118A/q 
models with a helium core composition and an envelope metallicity of Z = 0.001. OC means oxygen - carbon 
core composition. (Note that OC sequences for model more massive than IMq are indicated with the same 
notation). We also provide a column (T e ! ff ) with the effective temperature (in thousand K degrees) at which 
each sequence starts. 



3 THE GRIDS OF WHITE DWARF MODELS 

In this section we comment on the most important fea- 
tures of the grids. We computed evolutionary sequences 
with masses ranging from M = 0.1 Mq to M = 1.2 
M Q and metallicity Z = 0.001 and Z = 0. For models 
with M < 0.5 Mq we assume a pure helium core and for 
models with M > 0.45 Mq we assume a carbon - oxygen 
core with the chemical profile of Figure ^. We also vary 
the mass of the hydrogen envelope Mh within the range 
10 -12 < Mh/M < 4 x 10~ 3 and the mass of the helium 
layer Mue (in the case of carbon - oxygen white dwarfs) is 
taken to be Mu c /M — 10 -2 . In Tables 1 and 2 we summarize 
the main characteristics of all of our available evolutionary 
sequences. As stated earlier, models have been calculated 
in the framework of Canuto et al. (1996) theory of convec- 
tion. We have used OPAL opacity calculations supplemented 
with the Alexander & Ferguson (1994) molecular opacities 
(or with the Cox & Stewart 1970 tabulation for helium com- 
position) for low temperatures. In Table 1 and 2 we also give 
the T ft values at which each evolutionary sequence starts. 
In this regard, we emphasise once again that model ages cor- 
responding to the first stages of evolution are meaningless 
because they are strongly affected by the procedure we use 
to generate the initial models. At advanced ages, however, 
this is no longer relevant and age values are meaningful. The 
sequences have been evolved down to a stellar luminosity 
log L/Lq = -5. 

We begin by examining the time evolution of our mod- 
els. From the point of view of an age determination of the 
disc of our Galaxy from the observed space density of white 
dwarfs, the evolutionary times of white dwarfs as a function 
of mass represent obviously an important issue (see Wood 
1992 and references cited therein). In this regard, we feel it 
to be valuable to compare our cooling curves against those 
published by other authors. We elect those computed by 
Wood (1995) for models with pure oxygen core on the basis 
of OPAL radiative opacities. To this end, we have computed 
additional sequences for models with oxygen cores and with 
the same outer layer chemical stratification and metallicity 



as considered by Wood. The comparison is shown in Fig- 
ure ^ for 0.5 and 0.7 Mq models. Note the good agreement 
between the two set of calculations. At very low luminosities 
and especially for more massive models than those consid- 
ered in Figure ^, some divergency appears as a result in part 
of the different set of low - temperature radiative opacities 
employed. 

As we mentioned, white dwarf evolutionary times rep- 
resent a powerful tool for constraining the age of the disc 
of our Galaxy. Indeed, the existence of an abrupt falloff in 
the observed white dwarf luminosity function (see Liebert, 
Dahn & Monet 1988 and earlier references cited therein) 
has been interpreted in terms of a finite age of the disc of 
the Galaxy (D'Antona & Mazzitelli 1978). By fitting the 
observations with theoretical white dwarf luminosity func- 
tions, this interpretation was quantitatively explored by nu- 
merous investigators such as Winget et al. (1987), Iben & 
Laughlin (1989) and Wood (1992). Recently, Leggett et al. 
(1998) have substantially improved the determination of 
the observed luminosity function for cool white dwarfs. To 
compare with observations, we constructed integrated lumi- 
nosity functions from our evolutionary sequences. To this 
end, we follow the treatment presented in Iben & Laughlin 
(1989). Specifically, the space density of white dwarfs per 
unit of (, I = log(L/ Xq), is calculated from 

^ = -H>o I 4>{m) (^^-)m dm. (1) 

J mi 

Here, 0(m) is the Salpeter initial mass function of white 
dwarf progenitors with stellar mass m (which predicts that 
the created star distribution is proportional to 1/m 2 ' 35 ) and 
t coo i is the white dwarf cooling time at a given I, which is a 
function of the white dwarf mass M. mi and m a denote re- 
spectively the minimum and the maximum mass of the main 
sequence stars which contribute to the white dwarf space 
density at i. We take m s w 8 Mq (Wood 1992) and m* is 
obtained by solving the equation tMs(m) + t coo i(l,M) — td, 
where td is the assumed disc age. The pre - white dwarf evo- 
lutionary times tMs{m) are those of Iben & Laughlin (1989). 
As far as the initial(m) - final(M) mass relation is concerned, 
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Figure 4. Theoretical white dwarf luminosity functions (dotted 
lines) corresponding to our carbon - oxygen core, white dwarf 
models with a hydrogen envelope mass of Mjj/M= 10~ 4 and 
metallicity Z = 0. The curves, which correspond to assumed disc 
ages of 6 - 10 Gyr (at intervals of 1 Gyr), are compared to the 
observational data of Leggett et al. (1998) and they have been 
normalized to the observed white dwarf space density of 0.00339 
stars per cubic parsec. Note that the best fit to the dimmest white 
dwarfs observed correponds to a disc age of approximately 8 Gyr. 



Figure 3. Age - surface luminosity relation for 0.5 (upper panel) 
and 0.7 (lower panel) Mq white dwarf models with pure oxygen 
cores and hydrogen envelope mass of Mn/M= 10 -4 . We compare 
our results (solid lines) to those given by the Wood (1995) models 
(open circles) having the same stellar mass and chemical strati- 
fication as ours. The calculations are for a metallicity of Z = 0. 
At high luminosities, the discrepancy between the two set of cal- 
culations are due to the different procedure employed to generate 
the initial models. 

we use an exponential model: M = 0.40 e° 125 m (Wood 
1992). In deriving Equation (|I|), the star formation rate ip 
has been assumed to be constant. Finally, for each of the 
selected luminosity values, we calculate dt coo i/d£ at a given 
M by using linear interpolation between the dt coo i / dl values 
of the sequences that bracket M. The resulting luminosity 
functions for assumed disc ages of 6 - 10 Gyr are shown in 
Figure ^ (in Figure ^, the luminosity functions have been 
converted into intervals of bolometric magnitude Mtoi). It 
is worth mentioning that all of our theoretical curves have 
been normalized to the observed space density of 0.00339 
white dwarfs per cubic parsec (Leggett et al. 1998). The 
best fit to the coolest white dwarfs observed is obtained for 
assumed disc ages of ~ 8 Gyr, which is in agreement with 
the ages quoted by Leggett et al. on the basis of the Wood 
(1995) cooling sequences. 

It is worthy of comment that the age of helium white 
dwarf models depends on the mass of the hydrogen enve- 
lope. This is particularly true for models with very thick hy- 



drogen envelopes, for which hydrogen burning contributes 
substantially to the total luminosity, thus leading to a delay 
in cooling even down to very low T G ff . Unfortunately, the 
maximum mass of the hydrogen envelope resulting from bi- 
nary evolution is still an open question. Evidence favouring 
"thin " envelopes was presented by Iben & Tutukov (1986) 
from self - consistent binary evolutionary calculations. In- 
deed, these authors found that the hydrogen envelope re- 
maining at the top of their 0.3 Mq remnant after shell flash 
episodes is too small (Mh w 1.4 X 10 -4 Mq ) to sustain 
any further nuclear burning. Needless to say, a larger hydro- 
gen remnant would lead to longer evolutionary times. In the 
context of age determinations of millisecond pulsars with 
helium white dwarf companions, this fact is a clearly impor- 
tant one to be taken into account. To place this assertion 
on a more quantitative basis, we consider the pulsar PSR 
J1012+5307. The surface gravity and T e g for its low - mass 
helium white dwarf companion have been determined to be 
log g = 6.75 ± 0.07 and T cff = 8550 ± 25if (van Kerkwijk et 
al. 1996). At Trff « 85007^, we find that our 0.21 Mq helium 
white dwarf model with Mn/M = 2 x 10" 3 has log g = 6.82 
and age 0.44 Gyr, in good agreement with the Hansen & 
Phinney (1998b) predictions. The same fit to the T e g and 
gravity values would be achieved with a 0.213 Mq model 
with Mn/M ~ 6 x 10~ 3 . However, in this case hydrogen 
burning supplies 70 per cent of the surface luminosity and 
the model age becomes as high as 0.9 Gyr. 

Our main motivation for the publication of a detailed 
set of low - mass white dwarf models like the one presented 
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Table 3. Masses for low mass white dwarfs 



White Dwarf 


Source 


T eS (K) 




Ml Mr* (11 


Ml Mr* (2) 


Ml Mr* (3) 


0316 - 


- 345 


BSL 


14880 


7.61 


0.40 


0.414 


0.452 


0339 - 


- 523 


BSL 


13350 


7.47 


0.34 


0.354 


0.395 


0710 - 


- 741 


BSL 


18930 


7.45 


0.35 


0.377 


0.418 


0957 - 


666 


BRB 


27047 


7.285 


0.335 


0.371 


0.417 


1022 - 


- 050 


BRB 


14481 


7.483 


0.351 


0.364 


0.405 


1101 - 


- 364 


BSL 


13610 


7.38 


0.31 


0.325 


0.369 


1241 - 


h 010 


BSL 


24010 


7.22 


0.31 


0.342 


0.390 


1317 - 


- 453 


BSL 


14000 


7.43 


0.33 


0.343 


0.385 


1353 - 


- 409 


BSL 


23580 


7.54 


0.40 


0.427 


0.470 


1614 - 


h 136 


BSL 


22430 


7.34 


0.33 


0.361 


0.409 


1713 - 


h 332 


BSL 


22030 


7.40 


0.35 


0.376 


0.420 


1824 - 


- 040 


BRB 


14795 


7.608 


0.394 


0.411 


0.451 


2032 - 


- 188 


BSL 


18540 


7.48 


0.36 


0.385 


0.425 


2331 - 


r 290 


BSL 


27830 


7.50 


0.39 


0.431 


0.479 


2337 - 


760 


BRB 


14295 


7.507 


0.354 


0.372 


0.410 



T e ff , the surface gravity (g) and the stellar mass (1) of the objects were taken from Bergeron, Saffer & 
Liebert (BSL) (1992) and Bragaglia, Renzini & Bergeron (BRB) (1995). The next column (2) gives the 
stellar mass according to our helium - core models without a hydrogen envelope and the last column lists 
the stellar mass according to our helium - core models with a hydrogen envelope of Mh/M = 10 — 3 . 



here is, apart from the fact that little attention has been paid 
in the past to the study of this kind of objects, that both 
finite - temperature effects and hydrogen envelopes substan- 
tially modify the surface gravity values of zero - tempera- 
ture, helium - core - degenerate configurations. These fea- 
tures may turn out to be very important for the interpreta- 
tion of the recent and future observational data about low 
- mass white dwarfs. We think that the detailed low - mass 
models we computed here by employing the full scheme of 
stellar evolution theory may help such an endeavour. 

In the context of the foregoing paragraph we show in 
Figure |B| the surface gravity g (in cgs units) in terms of 
T e ff for some selected models. The effects of finite tempera- 
ture are clearly noticeable, particularly for less massive mod- 
els. As is well known, at a given T e g more massive mod- 
els are characterized by smaller radii. As T e g decreases the 
model radius (gravity) gradually becomes smaller (larger), 
ultimately reaching an almost constant value as expected for 
a strongly degenerate configuration, in which the mechanical 
structure is determined mainly by degenerate electron pres- 
sure. As a result, stellar parameters asymptotically reach 
constant values corresponding to zero temperature configu- 
rations. Note also the changes in the g values brought about 
by rather thick hydrogen envelopes. Another observation we 
can make from this figure is that, at low T e ff , convective 
mixing between hydrogen and helium layers increases the 
g values of models with thin hydrogen envelopes. In fact, 
convective mixing changes the outer layer composition from 
a hydrogen - dominated to a helium - dominated one, thus 
giving rise to denser outer layers. From then on, their sub- 
sequent evolution resembles that of a white dwarf model 
without a hydrogen envelope, as can be noted from Figure ^ 
(see also Figure pi), particularly for less massive models. 

To clarify better the role played by hydrogen envelopes 
in the g values of low - mass models, we show in Figure ^ 
the reduction in the g values of pure helium configurations 
resulting from adding hydrogen envelopes of different thick- 
ness. More precisely, we plot in terms of T e ff the quantity 
Ag/g = (g a — gn)/go for various stellar masses and hydro- 
gen envelopes (g and gn stand for the surface gravity of a 
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Figure 5. Surface gravities versus T c ff for selected white dwarf 
models with M/Mq = 0.15, 0.25, 0.35, 0.45, 0.60 and 0.80 and 
different hydrogen envelopes. For each stellar mass and from top 
to bottom the curves correspond to sequences having hydrogen 
envelopes with fractional mass of Mn/M= (no hydrogen en- 
velope), 10~ 6 and 10~ 4 , respectively. Note that in the case of 
less massive models, hydrogen envelopes appreciably reduce the 
surface gravity values of pure helium models. 



helium - core configuration of a given stellar mass without 
and with a hydrogen envelope, respectively). It is clear that 
thick hydrogen envelopes appreciably reduce the g values of 
pure helium models. At the low T e ff of 15,000K for instance, 
the g values of the 0.35 and 0.25 Mq pure helium models are 
reduced, respectively, by 20 and 30 per cent if a hydrogen 
envelope of Mh/M = 10~ 4 is added. Such values increase 
considerably at higher T e g . At T e s ~ 17,000 K, there is a 
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Ag/g 




T eff [x 10 3 K] 

Figure 6. (a): Ratio of the difference in g values between low 
- mass white dwarf models without (g ) and with hydrogen en- 
velopes to g versus T e g for models with M/Mq = 0.15 (solid 
lines) and 0.25 (dashed lines). For each stellar mass and from top 
to bottom the curves correspond to sequences with hydrogen en- 
velopes of fractional mass of M H /M= 4 X 10~ 4 , 10 -4 , 10~ 6 and 
10 -8 , respectively, (b): As (a) but for models with with M/Mq = 
0.35 (solid lines) and 0.45 (dashed lines). Note that thick hydrogen 
envelopes appreciably reduce the g values of pure helium models. 
Note also the effect of convective mixing at low T e fj on models 
with thin hydrogen envelopes. 

change in the slope of the curves stemming from the decrease 
in the radiative opacity values after helium recombination. 
This causes pure helium models to become denser and hence 
to have larger g values. 

Non - negligible differences in the structure and cool- 
ing of white dwarfs may also arise from the employment 
of different metallicities in the envelope, particularly at low 
luminosities where central temperature of models becomes 
strongly tied to the details of the outer layer chemical strat- 
ification (see Tassoul et al. 1990). This expectation is borne 
out by Figure 0, in which we compare the cooling times of 
helium white dwarf models for two extreme metallicities as- 
sumed in the envelope. When convection reaches the domain 
of degeneracy, the central temperature drops substantially 




2 4 6 8 10 12 

r 9 n 

Age [10 yrj 

Figure 7. Surface luminosity versus age relation for (from top to 
bottom) 0.40, 0.30 and 0.2 Mq helium white dwarf models with 
Mn/M= 10~ 4 and for metallicities Z = (dashed lines) and Z = 
0.001 (dotted lines). Note that at high luminosities, cooling is not 
affected by the assumed metallicity in the envelope. 

and the star has initially an excess of internal energy to be 
radiated, thus giving rise to a lengthening of evolutionary 
times during that stage of evolution. Because models with 
lower metallicities are characterized by deeper convection 
zones (see Figure [l]) , this effect occurs at higher luminosities 
in such models and this explains their greater ages as com- 
pared with high - metallicity models. Eventually, at very low 
luminosities, more transparent models evolve more rapidly, 
as expected. We have also analysed the effect of metallicity 
on surface gravity for helium models (see Figure ^|) and we 
found that surface gravity is almost insensitive to a specific 
choice of metallicity in the envelope. 

Lastly, we have applied our evolutionary models with 
helium cores to derive stellar masses of some selected low - 
mass white dwarfs. To this end, we have picked out low sur- 
face gravity white dwarfs from the sample of white dwarfs 
analysed by Bergeron, Saffer & Liebert (1992) and Bra- 
gaglia, Renzini & Bergeron (1995), and we list the results in 
Table 3. The above cited authors estimated the white dwarf 
masses from evolutionary models with pure carbon - core 
composition. However, using our evolutionary models with 
helium cores and no hydrogen envelope we find the mass val- 
ues to be appreciably underestimated, particularly at high 
temperatures. It is worth mentioning that the objects listed 
in Table 3 are white dwarfs that most likely have a hydrogen 
envelope. This being the case, the stellar mass should be es- 
timated from evolutionary models with hydrogen envelopes. 
As shown in Table 3, there is an appreciable difference in 
the white dwarf mass when stellar masses are derived from 
models with thick hydrogen envelopes (Mh/M « 10 3 ). 

In view of the preceding considerations, we conclude 
that detailed models of low - mass white dwarfs such as pre- 
sented in this study should be carefully taken into account, 
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Figure 8. Surface gravities versus T e ff for (from top to bottom) 
0.30, 0.25 and 0.20 Mq helium white dwarf models with M H /M= 
and for metallicities Z = (solid lines) and Z = 0.001 
(dashed lines). 



should the mass of a white dwarf be measured by apply- 
ing the surface gravity - T c s relation. This is particularly 
true regarding the possibility of constraining the equation 
of state at neutron star densities as inferred from observa- 
tions of low - mass white dwarf companions to millisecond 
pulsars, such as those studied by van Kerkwijk et al. (1996). 

Complete tables containing the results of our cal- 
culations are available at the World Wide Web site 

Additional evo- 



http : // www. fcaglp . unlp .edu .ar / ~ althaus / 



lutionary sequences are obtained upon request to the au- 
thors at their e-mail address. Features such as surface lumi- 
nosity, T ff , central density and temperature, surface grav- 
ity, stellar radius, age and hydrogen surface abundance are 
listed in the tables. 
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